Preprint typeset in JHEP style - PAPER VERSION 



October 2012 
FPAUO-12/14 



Af = 8 SYM vs. Af = 6 Chern-Simons: 
Four-point amplitudes at two-loops 



Marco S. Bianchi and Matias Leoni' 

* Departamento de Fisica, Universidad de Oviedo, 33007, Oviedo, Spain 
Pi 

^Physics Department, FCEyN-UBA & IFIBA-CONICET Ciudad Universitaria, 



Pabellon I, 14-28, Buenos Aires, Argentina 

E-mail: bianchimarcoOuniovi . es , leoniSdf . uba . ar 

Abstract: We compute the four-point amplitude in 3d J\f — 8 SYM at two loops, 
by solving the three dimensional scalar doublebox in dimensional regularization. We 
compare it to the same result in the ABJM theory, to which maximal SYM should 
flow in the infrared at strong coupling. After proper rescalings, we find that the two 
amplitudes coincide in the Regge limit. 
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1. Introduction. 

Amazing progress in the computation of scattering amplitudes and in a deeper under- 
standing of their properties has been attained in the last years. Most of these advances 
have been achieved in the realm of M = 4 SYM theory in four dimensions, however, in 
order to gain a broader insight on the subject, it is worthwhile exploring amplitudes in 
other contexts. Since the realization of the ABJM model [1, 2] in three dimensions, an 
exciting new environment where one can study amplitudes was opened. Many results 
are now already available for Af = 6 Chern-Simons matter theories: after the first in- 
vestigation of [3], the structure of tree level amplitudes was analyzed [4, 5], discovering 
an underlying dual superconformal symmetry [6], which combined with the original 
superconformal invariance of the model, gives rise to a Yangian [7]. This symmetry has 
been proven to extend to all tree level amplitudes and loop integrands [8], by virtue of 
a three dimensional version of recursive relations [9, 10], and was exploited to propose 
an orthogonal Grassmannian integral formula [11], in close relationship with the M = 4 
one [12]. 

While this appears to be a symmetry at weak coupling, a neat dual interpretation, 
like the one of [13, 14] in four dimensions, is still lacking [15]— [21]. 
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At higher order, amplitudes have been computed at one loop through unitarity 
[22]-[24] and directly through Feynman supergraphs [25]. Such amplitudes are IR 
finite and show a three dimensional version of the holomorphic anomaly [26]. 

At two loops the four-point amplitude was computed [27]-[29]. It exhibits dual 
conformal invariance [30], duality with the bosonic four-cusped Wilson loop [31] com- 
puted in [32] and there are hints that it may enjoy exponentiation, similarly to the BDS 
four dimensional formula [33]. The two-loop four-point amplitude in ABJM strikingly 
resembles its one-loop analogue in four dimensional M = 4 SYM, and a relation con- 
necting the two, to all orders in the dimensional regularization parameter expansion 
was determined in [34], which might potentially extend to all loops, if the amplitudes 
do really exponentiate. 

Furthermore the computation of the six-point amplitude at two-loops has been 
recently carried out [35]. 

Apart from Chern-Simons theory, amplitudes have been studied in three dimen- 
sions for supersymmetric Yang-Mills models. The Yang-Mills action being super- 
renormalizable and having a dimensionful coupling, prevents the theory from being 
superconformal, and therefore from having an AdS/CFT dual. Hence one may ex- 
pect that the nice properties of ABJM amplitudes may not be shared by SYM ones. 
Still, these theories may reveal interesting aspects, since they are expected to flow to 
Chern-Simons matter theories in the infrared. 

In particular maximal supersymmetric Yang-Mills theory with gauge group U(N), 
which describes the low-energy worldvolume theory of a stack of N D2-branes, is 
believed to have ABJM U(N)i x U(N)-i as a conformal fixed point of its RG flow in 
the deep infrared. In this phase the analysis is hampered by the inherently strongly 
coupled regime of the two theories, nevertheless evidence of their equivalence has been 
pointed out in [36]-[38], matching their partition functions and superconformal indices. 

It is not clear how this may reflect on the structure of scattering amplitudes, how- 
ever [39] showed how the 5*0(8) R-symmetry expected for Af = 8 enhanced ABJM 
with level one can be recovered at four points in maximally SYM scattering, and such 
an analysis was extended in [40]. 

By analyzing tree level amplitudes of Af = 8 SYM in on-shell superspace, and by 
means of recursion relations, [41] proved that they enjoy dual conformal covariance. 
This also applies to the integrands of loop amplitudes by means of unitarity. Further- 
more, that paper developed a prescription to go from four dimensional SYM amplitudes 
to three dimensional ones by means of dimensional reduction, and used it to show that 
Af = 8 SYM amplitudes in three dimensions have helicity structure and that beyond 
four-points, the amplitudes should have SO (7) R-symmetry. 
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In the same paper, the authors also argued that in maximal three dimensional 
SYM, the loop integrands can be obtained from four dimensional M = 4 SYM and 
subsequently showed that at one loop, all MHV amplitudes vanish and all non-MHV 
amplitudes are finite. This is at least consistent with what occurs in ABJM. The 
authors in [41] also point out that at two loops the four-point amplitudes of ABJM 
and Af = 8 SYM might show some similarity, motivated by their conjectured relation 
and by the fact that they share dual conformal covariance. 

In this short note we explore this possibility by explicitly computing the two-loop 
ratio of the four-point amplitude to its tree level counterpart for M = 8 SYM in three 
dimensions. 

Following dimensional reduction, this task may be accomplished by computing the 
three dimensional scalar doublebox. As expected, this integral suffers from infrared 
divergences which we regulate through dimensional regularization. This computation 
is described in section 2, where we also give the final expression for the amplitude. 

In section 3 we comment on the outcome. Given the dimensionful nature of the 
coupling constant, we expect the ratio to have mass dimension —2 at two-loops, which 
is indeed the case, due to an overall (s~ x + t^ 1 ) factor. The amplitude presents infrared 
singularities, in the form of e~ 2 poles, whereas subleading divergences may be absorbed 
into a proper rescaling of the mass parameter of dimensional regularization. The finite 
part is made of a squared logarithm of the ratio of Mandelstam invariants, a subleading 
logarithm multiplied by an s -H- 1 antisymmetric combination, and a constant piece. In 
contrast to ABJM, this amplitude does not exhibit uniform transcendentality. 

Since the theory does not even possess conformal invariance one suspects that a 
duality with light-like Wilson loops should not occur. We verify that this suspicion is 
correct by computing the first order correction to the four cusped Wilson-loop, which 
is non-vanishing in contrast to the result for the amplitude of [41]. At two loops we 
also check that not even the UV divergences of the cusp in the light-like Wilson loop 
match the IR singularities of the amplitude. Actually at two loops the Wilson loop only 
shows e _1 UV poles, as a result of milder short distance singularity, since the theory is 
superrenormalizable. 

In the last section we compare the SYM amplitude with the ABJM one. Upon 
redefining the relative coupling constants in such a way that the effective YM param- 
eter is dimensionless, and that the coefficients in front of the infrared divergent piece 
coincide, the two amplitudes exhibit some partial resemblance, spoiled by the presence 
of the term proportional to log s/t, absent in the ABJM case. Quite interestingly, the 
maximal transcendentality part of the constant coincides with that of ABJM. Taking 
the Regge limit s/t — > 0, the leading logarithm approximation reproduces exactly the 
ABJM case, in the sense that the log 2 s/t coefficient is precisely the same. For ABJM, 
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thanks to the duality with the Wilson loop, an anomalous conformal Ward identity 
[42, 43] fixes this coefficient to be one half that multiplying the — e~ 2 poles. 



2. Computation of the amplitude: the three dimensional scalar 
doublebox. 

We consider Af = 8 SYM theory in three dimensions, with unitary gauge group U(N), 
whose rank we take large in the planar limit N ^> 1, and coupling constant gyM, having 
mass dimension 1/2. 

At weak coupling, and in the planar limit, we expand the color ordered four gluon 
amplitude in the power series 

where e = is the dimensional regularization parameter and /x a mass scale. 

We want to compute the two-loop correction to M4. Following the derivation 
in [41], this amounts to borrowing the four dimensional result for the integrand and 
shifting the dimension of the Feynman integrals to d = 3 — 2e. This is motivated 
by the observation that the integrands of three dimensional Af = 8 SYM obey the 
same transformation properties of those of Af = 4 SYM in four dimensions, under dual 
inversion, and are thus the same. 

Therefore we compute the three dimensional scalar massless doublebox 1 in Fig. 1. 



d 3 -^kd 3 ~ 2 n 

1 ' " ! k 2 (k - p 2 ) 2 (k + Pl ) 2 (k + l + Pl + Pi ) 2 (l + - p 3 ) 2 l 2 



X (2) = _ e 2 7B£7r -3+2 £ 



-3+2e / 

* J 



(2.2) 

Note that this integral is not dual conformally invariant, a property which is prevented 
because we have lowered the dimension of the integration measure to three. However it 
inherits dual conformal covariance from four dimensions. Calling s/t = x, its Mellin- 
Barnes representation reads 



1 { £ = -T^JM ( 2 -3) 



1 We follow the conventions of [33] for the normalization of the integral. 
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Figure 1: The scalar doublebox. 



where 

1 p+ioo 4 j 

r (—^ - -23 - 24) r + Z4) r + z 4 ) r (z l + z 2 + z 3 + z 4 + 1) x 

r (-e - Zx -z 2 -\)Y (-e + z 2 - \) V (-e - z x - z 3 - §) T (-e + z 3 - §) X 

r( e + ^-^ 4 + |)r(£ + z4 + |) s-' 1 f , 

r (22 + z 4 + 1) r (z 3 + z 4 + 1) r (-2e + z x - z 4 ) 

The complex contour of integration in this four-fold Mellin-Barnes integral takes all 
variables in a straight line from —ioo to zoo with possible indentations, if necessary 
in order to leave the whole series of poles coming from Gamma functions of the form 

+ z) to the left of the contour and the whole series of poles coming from Gamma 
functions of the form T(... — z) to the right. Note that the contour of integration is well 
defined thanks to the presence of the e regulator inside the Gamma functions since it 
separates left and right poles which would otherwise collide. 

Since there is a F(— 1 — 2e) in the denominator of (2.4), the expression could seem 
of order e at first sight; this is a common feature in the Mellin-Barnes representation 
of Feynman integrals. Actually, due the phenomenon of colliding poles in the e — > 
limit, the integral itself has poles in e producing a non vanishing result. The way of 
proceeding is to analyze case by case where this situation is produced and deforming 
the contour by picking up residues so as to avoid potentially colliding poles in the e — > 
expansion. Thanks to this method, the Mellin-Barnes integral becomes effectively at 
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most two-fold, since we discard (9(e) terms. This procedure is explained in full detail 
in the literature (see for instance [44]) and has been successfully applied to many 
four dimensional examples [45]. In this case we employed the package MB.m [46], to 
automatically perform the e expansion and get a list of Mellin-Barnes integrals whose 
integrands are independent of the regularization parameter. From this analysis, one 
can infer the functional form of the result up to order e° 



T (2) vr e 2 ^ 

-La 



(_ s )4+2e 



(ai x 2 + a 2 x) + i ( (a 3 x 2 + a 4 x) log(x) + a 5 x 2 + a 6 x) + 



+ [a 7 x 2 + a 8 x) log(x) + a 9 x 2 + a w x 



+ 0(e) (2.5) 



where again x = s/t is the ratio between the Mandelstam kinematic variables and 
the n factor has been collected for future convenience. It is straightforward to derive 
this functional dependence on x, since due to the procedure of picking up poles in e 
described above, there is no explicit dependence on x in the Mellin-Barnes integrals, 
which therefore only contribute to give numerical coefficients, but not functions of 
the kinematic variables. The nontrivial logarithmic dependence of the result is indeed 
produced just by the e expansion. 

After the expansion is performed, we obtain a total of 108 one-fold and two-fold 
integrals (a list may be found in the attached Mathematica notebook) needed to get 
the values of the coefficients a«. This computation is a little bit cumbersome, due to the 
large amount of integrals, but in every case, after appropriate manipulations, they can 
all be reduced to a form where some corollary of the Barnes' first and second lemmas 
may be applied. Almost all of the 108 integrals could be solved with the extensive list 
of Barnes' lemmas corollaries detailed in [44]. Those which were not present in that 
list, are explained in Appendix B. 

Following this procedure we analytically determined the value of the coefficients to 

be 

ai — 6 , a 2 = 4 , a 3 = 6 , a 4 = 4, 
a 5 = 20 - 12 7s + 12 log 2 , a 6 = 27 - 8-f E + 8 log 2, 
a 7 = 28 - 12 7i? , a 8 = 28 -8 7 £, 
a 9 = -28 - 407b + 12 7 1 - 8tt 2 - 24 log 2 2 + 80 log 2 - 24 7s log 2, 

167T 2 

a 10 = 46 - 54 7e + 8 7 | - 16 log 2 2 + 52 log 2 - 16 7e log 2 (2.6) 

Plugging these into (2.5) one gets the three dimensional scalar doublebox (in the s- 
channel) up to 0(e). 
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In order to compute the whole four-point amplitude, we use the four dimensional 
result [47] and plug the three dimensional integrals we compute 



Mf (.si;.; -VtT 



'(2) 



At J 



(2.7) 



where the latter integral can be obtained from the former by replacing s <H> t. After 
some algebra and suitably redefining the mass scale /i as 



V = V exp ( — + log2 



47 



(2.8) 



we obtain the form of the two-loop 3d M = 8 four-point amplitude, which constitutes 
the main result of the paper 



Mf = 5tt 



s + 1 

St 



-8/fJ, 



/\-2e 



:-t//iO 



-2c 



+ ^log 2 (^)+4C 2 + 31og 2 2 



1 /, n 7\ s- 1 , /s\ 19 , n /43 



(2.9) 



3. Properties of the result. 

We comment some relevant properties of the result. First we note that since M. is the 
coefficient in the perturbative series of powers in the coupling constant g\ M , which has 
mass dimension 1, the two-loop ratio is dimensionful as well. The overall scale (- + 7) 
presents single poles in the collinear regimes s — > and t — > 0, in contrast with ABJM 
[27, 28]. 

Second we verify that infrared divergences appear at two loops. These are captured 
by the 1/e 2 poles, whereas the subleading pole may be non-trivially reabsorbed by the 
scheme change (2.8). As expected from the planar limit, these divergences come from 
an s channel and a t channel contributions. After reabsorbing the 1/e pole the finite 
piece of the amplitude becomes independent of the \i regularization scale. 

Comparison to Wilson loops. As the theory obviously lacks both ordinary and 
dual conformal invariance, one should not expect a WL/amplitude duality to work. 
To ascertain this, we compute a four cusp light-like Wilson loop at one loop. The 
necessary tools for this are given in Appendix A. In order to get a real result we 
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restrict to a light-like polygonal contour with space-like non-vanishing invariants (the 
diagonals), and set s = xj 3 < and t = x| 4 < 0. We also remark that at one loop 
no UV divergences arise, so that the calculation may be safely performed in d — 3 
dimensions. Finally the expectation value of this Wilson loop reads 

= , _ ArcTanh / V^ + V^T) \ + a(A ) (3 ,, 

7T a/ — S — t \ y/st — S — t I 

where CV = T a T a is the quadratic Casimir of U(N), which is N/2 in the fundamental 
representation. Comparing it with the amplitude result, which vanishes at this order, 
we verify that the duality already fails in the simplest case. 

Not even the UV divergent piece of the light-like Wilson loop, which appears first at 
two loops, resembles the IR divergences of the amplitude. We checked this by analyzing 
planar corrections to a Wilson light-like line cusp at second order in perturbation 
theory, from which the divergences of the polygonal Wilson loop arise. The computation 
parallels the four dimensional one and we won't go through it, but just state the result. 
The relevant tools for it to be performed may be found in Appendix A. We verified that 
by virtue of the milder UV behavior of the gluon x-space propagator (A. 2) in three 
dimensions, where YM theory is superrenormalizable, the rainbow diagram 2 as well as 
those involving the cubic gluon interaction are finite. The only source of divergence 
comes from the gluon self energy insertion. The one-loop corrected propagator (A. 7) is 
indeed finite in momentum space. However, on dimensional grounds, its dependence on 
p 2 is through (— p 2 )~ 3 ^ 2 ~ e . In dimensional regularization, its Fourier transform develops 
a single pole in e, which accounts for the UV divergence. In fact, once the corrected 
propagator is inserted in the Wilson loop, no new singularities appear and the divergent 
piece of the Wilson loop reads 

WOk- -2*^4 <-»-«) (3.2) 
where s and t stand again for the diagonals of the WL contour. 

Transcendentality. We observe that the result does not respect maximal transcen- 
dentality, in contrast with ABJM. Once the integral is normalized as in (2.2), all ^ E 
factors disappear, although leaving coefficients with mixed transcendentality, such as 
the one multiplying the log s/t piece, or the constant. Still, we point out that the 
maximal transcendentality piece of the constant coincides with that of ABJM. 

2 The one obtained by expanding the Wilson loop exponential to fourth order and Wick contracting 
two pairs of gluons lying on opposite edges of the cusp in a planar way. 
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4. Comparison to ABJM: the Regge limit. 

We want to make a comparison of the M = 8 SYM two-loop four-point ratio (2.9) 
with the same result for ABJM, which we recall 

Mf — = - <Z£p - Hp + \ V* g) + 4 & + 3 W 2 (4.1) 

The form of this amplitude coincide, up to rescalings and up to the constant, with the 
one-loop amplitude of M = 4 SYM in four dimensions 

M ^ = J-fl-tfl + i W (Vj +i(2 (4.2) 

Comparing (2.9) with (4.1), we find the following relation 

n-^Mf )M = 8 = 5Mf )ABJM +(\o g 2-l) ^ log"- + C + 0(e) (4.3) 

Curiously enough, written in this form, the part of the constant of the M — 8 SYM 
result with maximal transcendentality coincides with the ABJM one and the difference 
are terms with lower transcendentality, namely 

19 43 2 
C = — -log2 + — (4.4) 
2 s 80 V ; 

The relevant part of the amplitude which differs from the ABJM result is the extra 
piece proportional to the logarithm of the ratio of the Mandelstam variables. This 
suggests that we might compare the M = 8 result with that of ABJM in the Regge 
limit t/s <C 1, where the leading logarithm approximation is pursued. To make the 
comparison sensible, we need to deal with objects of the same dimension. While the 
perturbative expansion of ABJM amplitudes is carried out in terms of the dimensionless 
't Hooft parameter \ C s = N/k, the Yang-Mills effective coupling constant in the 
planar limit is \ym = 9ymN, which in three dimensions is dimensionful, having mass 
dimension one. While the Chern-Simons coupling may be tuned to be small adjusting 
the rank and the level in such a way that N <^ k, in SYM the coupling constant 
has to be compared to another energy scale. In the four-point amplitude there are 
two characteristic energies of the process associated to t and s Mandelstam variables. 
Therefore perturbation theory is valid whenever g YM <^ y/s, \/t. We may further fix t 
to be the energy scale we are interested in, and study the small t/s regime. In such a 
limit we can rescale the Yang-Mills coupling g YM = y/i, defining a dimensionless 
effective coupling g 2 e ^ = gyul^-, which we take to be small g e ff 2 <C 1. In this way we 
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establish the hierarchy of energy scales we are interested in, which is g\ M <C yt <C \fs. 
In addition, to make the comparison more direct, we can focus on the finite pieces of 
the SYM and ABJM amplitudes, rescaling the latter by the 57r factor, such that the 
normalization of the infrared divergences coincides between the two. Having done this, 
and in terms of the dimensionless coupling, the M = 8 amplitude becomes 



finite 



;i + y) 



- log 2 y + 4C 2 + 31og 2 2 



1 



1 



^(log2-|) ^—^ logy- flog2+ (f ) 2 



(4.5) 



where y = l/x = t/s and A 



eff 



9lff N - 




Figure 2: Log-Log plot of M = 8 and ABJM amplitudes as functions of y = t/s. 
In the Regge limit y 1 this becomes 



Mf )N=8 



finite 



^Ilog 2 y + 0(logy) 



while ABJM, in the same limit becomes 

M 



(2) ABJM I y«l 1 , 2 , /ntl\ 

->■ g lo s y + c> ( 1 ) 



I finite 



(4.6) 



(4.7) 



Ergo, the finite piece of both amplitudes coincide in the Regge limit if we identify their 
dimensionless 't Hooft couplings A^ 5 4-> A 2 jj. In Fig. 2 we plot equation (4.5) and 
the finite piece of (4.1) as functions of y = t/s in logarithmic scale to illustrate the 
coincident Regge limit for small t/s. 

Amplitudes in = 6 CSM theory exhibit dual conformal invariance and, in the 
four point case, duality to the light-like four-polygonal Wilson loop. This has dramatic 
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consequences on the dependence of the amplitude on the kinematic invariants, which 
has to satisfy an anomalous conformal Ward identity. In the four-point case, at two 
loops, this fixes the finite piece of the amplitude to have a log 2 s/t form and its coefficient 
to be half that governing the IR singularities. It is thus very pleasing that in Regge 
limit the leading term of the finite piece of the MSYM amplitude is exactly log 2 s/t, 
and its coefficient of the form dictated by the anomalous conformal Ward identity. 
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A. Tools for the WL computation. 

We work in three dimensional Minkowski spacetime with metric g^ v = diag (1,-1,-1). 
We use dimensional regularization of momentum integrals d = 3 — 2e, and introduce a 
mass scale ji to keep the action dimensionless. The Wilson line operator is defined as 



W[C] = —TrVe ig f cdxtlA ^ 



(A.l) 



The x-space propagator for the gluon in 3c? SYM, in dimensional regularization and in 
the Feynman gauge, reads 



(A ll (z)M0)) = GV>(z) = - "* fJ "» 



-x 



2NV2-e 



(A.2) 



The momentum space one-loop corrected gauge propagator reads 



(A,( P ) m-p)) {1) = g$(p) = ( 4 ) ng(p) ( 



(A.3) 



where the gluon polarization operator evaluates, in dimensional regularization, 



; r(l- e )r(|-e)r(| + e ) 



((7-6e)JV-(2-4e)n/-n a ) 



(47r) 3 / 2 - e /i- 2e r(3-2e) 
For Af = 8 SYM, where n f = AN and n s = 7N, 

^ W ^ ^ (4vr)3/2-r(3-2e) 

and Fourier transforming to x-space with 

d d p e~ ipx i r (f - a) 



P Qjiv - VyPv 
(_ p 2)l/2+e 

(A.4) 



^2 



(27r) d (-p 2 ) a 4 a 7r rf / 2 r(a) 



(_ a .2 ) a-d/2 



(A.5) 



(A.6) 



this becomes 



2c 



d 3 2e p 



3-2e 



r (~ 2£ ) (_ x 2 )2e _ I r(-l-2e) " 

r(| + e ) 1 j ^ 4 r(| + e) ^ l j 



(A.7) 
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B. New corollaries. 



To complete the computation we also had to derive two corollaries of Barnes' first 
lemma, to solve the integrals (in the first Ai + A 3 7^ is assumed) 



-1 rioo 

— dzT(z + Ai) r (z + A 2 ) r* (-z - A 2 - 1) r (A 3 - z) ^ (0) (-z -X2-I), 
2m J_ ioo 

— / dzT (z - Ai + 1) r* (z - A 2 - 1) T* (A 2 -z)Y (Ax - 2) V (0) (A 2 - 2) 

2tt2 ./_ i00 

These may be found from lemmas (D.12) and (D.14) of [44]. To do this one may first 
insert a fictitious regulator to avoid coincident left and right poles. Having done this 
one takes a residue to shift the contour in such a way to remove one asterisk from the 
proper Gamma function and bring the integral to a form where the aforementioned 
corollaries of first Barnes' lemma may be applied. Since the original integral was well- 
defined, the poles in the regulator disappear in the sum of the two pieces and the limit 
where it approaches zero can be safely taken, to get the final result. We thus obtain 

1 rioo 

— dzT(z + X 1 )T(z + A 2 ) r* (-z - A 2 - 1) T (A 3 - z) ^ (0) (-z - A 2 - 1) = 
2m J_ ioo 

12r(A 1 + A 3 -l) P (Al " A2) 1 (Al + A3 " 1} 1 (A2 + A3) X 
(-12 7 ^ (0) (Ax - A 2 ) + 12 7 ^ (0) (A 2 + A 3 ) + 7T 2 + 6 7 2 ) + 

+ r (Ax - a 2 - 1) (r (Ax + a 3 - 1) r (a 2 + a 3 + 1) (-12(7 - i)v (0) (x 1 - a 2 - 1) + 

+12(7 - l)^ (0) (A 2 + A 3 + 1) + 7t 2 + 6 7 2 - 12 7 + 12) - r (Ax + A 3 ) r (A 2 + A 3 ) x 
(6^ (0) (Ax + A 3 - 1) 2 + 12( 7 - 1)V (0) (Ax + A 3 - 1) - 6^ (0) (A 2 + A 3 ) 2 + 

— 12V (0) (Ax - A 2 - 1) (V> (0) (Ax + A 3 - 1) - ^ (0) (A 2 + A 3 ) + 7 - 1) + 

-6^ (1) (Ax + A 3 - 1) - 6V> (1) (A 2 + A 3 ) + 7T 2 + 6 7 2 - 12 7 + 12))] (B.l) 

— / dzT(z-X 1 + l)T*(z-X 2 -l)T*(X 2 -z)T(\ 1 -z)^(\ 2 -z) = 
2m J_ ioo 

^[r (Ax-A 2 )r(-Ax + A 2 + l)x 

(12(7 - 1)V> (0) (Ax - A 2 ) - 12(7 - 1)^ (0) (-Ax + A 2 + 1) + 7t 2 + 6 7 2 - 12 7 + 12) + 

- T (Ax - A 2 - 1) (I2r (-Ax + A 2 + 1) (^ (0) (-Ax + A 2 + 1) + l) + 

+ T (-Ax + A 2 + 2) (6^ (0) (Ax - A 2 - 1) 2 - 12^ (0) (-Ax + A 2 + 2) ^ (0) (Ax - A 2 - 1) + 
+ 6^ (0) (-Ax + A 2 + 2) 2 + 6^ (1) (Ax - A 2 - 1) + 6^ (1) (-Ax + A 2 + 2) + 
-tt 2 - 67 2 + 127 - 12))] (B.2) 
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